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1. Introduction 

The problem of contraction of an air bubble in a Hele-Shaw cell filled with 
a Newtonian fluid under the influence of suction of air from the bubble has 
been intensively studied by many physicists and mathematicians for more 
than 25 years, see |GV ] and references therein. In particular, in |EElj . the 
authors suggested an analytic theory which allows one to give a complete 
description of the asymptotics of the bubble shape as its area goes to zero, 
and, in particular, to find the point of its contraction. In a number of 
applications of the contraction problem, for instance, in the theory of gas 
recovery, the following question, only briefly discussed in |EElj . is important: 
will the air bubble fall apart during contraction, or will it remain connected 
until all the air has been extracted? In this paper, we study this question 
in detail. In particular, we propose some sufficient conditions of breakup 
of the bubble, and ways to find the contraction points of its parts. In the 
theory of gas recovery, these points are interpreted as the optimal positions 
of the gas-producing wells. 

We note that all the results of this paper (with the exception of explicit 
solutions) extend to flows in a curved Hele-Shaw cell, along the lines of 
[EE21IEV] . 

The structure of the paper is as follows. In Section 2, we describe the 
mathematical model of the contraction problem, and recall some of the 
results from [EEl] on the connection between the dynamics of the bubble 
and its gravity potential. We also study the asymptotic shapes of bubbles 
contracting to a degenerate critical point of the potential. In Section 3, 
we define the points of partial contraction, which are contraction points of 
the bubbles which appear as a result of breaking of the initial bubble, and 
contract before the full contraction occurs, and extend to them the results 
about points of complete contraction from |EElj . In Section 4, we give a 
sufficient condition of breakup of a symmetric bubble; this is an extension 
of a result from [EEl ]. In Section 5, we consider the process of regulated 
contraction, which is simultaneous contraction of two or more bubbles under 
prescribed rates of extraction from each bubble. In Section 6, we consider 
the case of a bubble which breaks up into two bubbles under contraction, 
and discuss the question whether one can use regulated contraction to make 
these two bubbles contract simultaneously (development of singlularities in 



solutions may be a problem). A strategy of extraction which allows one 
to do so is called a synchronizing strategy, and we study such strategies 
in some detail. In particular, in Section 7 we study the asymptotics of 
contraction under a synchronizing strategy, and show that, like in the case 
of a single bubble, the two bubbles contract at critical points of the potential 
(generically, nondegenerate local minima), and the bubble shapes generically 
tend to ellipses, whose axes are determined by the eigenvalues of the Hessian 
of the potential at the minima. In Section 8, we characterize domains that 
are on the boundary between those that admit a synchronizing strategy and 
those that don't: we show that the potential of such a domain should have 
a degenerate critical point of the potential; we also study the asymptotics of 
contraction to such a point. In Section 9 we characterize domains that are on 
the boundary between those that break up and those that don't; generically 
they develop an instantaneous 5/2-cusp in the process of contraction. In 
Section 10 we correct some computational errors in the previous publications 
[EETl lEV]. 

Acknowledgments. The work of the second author was partially sup- 
ported by the NSF grant DMS-0504847. He thanks B. Gustafsson and D. 
Jerison for useful discussions. 

Note. Sadly, the first author of this paper, Vladimir M. Entov, passed 
away on April 10, 2008. 

2. The mathematical model and its main properties 

2.1. The model. Let us recall the formulation of the problem and the main 
results from the paper [EElj . Consider contraction of an air bubble in an 
unbounded Hele-Shaw cell, which is filled with a Newtonian fluid, under 
suction of air from the bubble. Let B{t) be the air domain at a time t. We 
assume that it is connected and bounded, with (say) a smooth boundary. 
Its law of evolution in time is as follows. In the fluid domain B{tY (the 
complement of B(t)), which we assume to be connected, there is a potential 
vector field of fluid velocities, v{x,y,t) = V^{x,y,t). The potential $ is 
determined at any time i as a solution of the boundary value problem 
(1^) 

l^+l^ = 0, {x, y) G B{tY- ^QB(t) = 0; $(x, y) = log(r)+0(l), r ^ 

where r = (x^ + y'^Y^'^, and (7 > is the rate of suction. The velocity of the 
boundary is then equal to the velocity of the fluid particles on the boundary: 

The contraction problem is to find the family of domains B{t) for a given 
initial shape of the bubble B{0) and given rate q of suction. It can be shown 
that for every initial domain with a smooth boundary, the solution of this 
problem exists and is unique on some interval of time [0, r), r > 0. It has 
an obvious monotonicity property B{ti) C -6(^2) for ti > t2- 



In a similar way one can define contraction of several bubbles. In this 
case, the domain B{t) is a union of finitely many disjoint simply connected 
domains (bubbles). 

2.2. Weak solutions. Equations ()l|2p determine contraction as long as it 
reduces to continuous deformation of the boundaries of the bubbles. How- 
ever, in the process of contraction, the boundary B(t) may undergo topo- 
logical transformations. For instance, parts of the boundary can collide 



In this case, the above definition of the contraction does not apply as it 
is, and needs clarification. 

Namely, let S be the area of B{0), and t* = S/q he the time of complete 
extraction of the air. Let us call a family of bounded domains B{t), t S 
[0,t*), a weak solution of the contraction problem if for every t E [0,t*), 
B(t) is a disjoint union of finitely many simply connected domains, so that 

(i) B{ti) C B{t2) for h > t2; 

(ii) the area of B{t) is 5 — qt; and 

(iii) there exists a closed set of times T C [0,t*] (of topological transfor- 
mations) containing 0, t* such that for any t £ T, a small enough interval 
{t, t + e) does not intersect T, and on the intervals of time {ti,T2) not inter- 
secting with T the domain B{t) is a "classical" solution of the contraction 
problem, i.e. is defined by equations (jl|2p . 

In other words, a weak solution is "glued" from usual (classical) solutions 
at points t £ T where the domain B{t) undergoes topological transfor- 
mations. It is natural to assume that it describes the actual process of 
contraction in the case when the breakup of the bubble does occur. 

It is known (see [GV]) that the contraction problem has a unique weak 
solution. An example of a weak solution is given on Fig.l; from this example 
one can see two types of topological transformations that can occur during 
contraction: 

1) breakup of a bubble into two pieces, and 

2) disappearance of a bubble. 

As a result of these transformations, the number of bubbles changes in 
the process of contraction. 

Remark 2.1. It can be shown that the set T is finite. However, a detailed 
proof of this would be long and we do not give it here. 

2.3. The gravity potential. Let us define the gravity potential of a bounded 
domain B to be the function 



where ^ = ^ + irj, z = x + iy. This function satisfies the Poisson differential 



(Fig.l). 



(3) 
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where Xsi^jV) is the characteristic function of B, and S is the area of B. 
It is shown in [EElj that this function is closely related to the contraction 
problem. Namely, we have the following theorem. 

Theorem 2.2. Let B{t) be a weak solution to the contraction problem. Then 
(i) The gravity potential inside B{t) changes by a constant under contrac- 
tion: 

(5) nB(o)(e,r/)-nB(,)(e,r/) = C(t), i^,rj) e Bit), 
(a) Let us set ^(x,y,t) = if {x,y) £ B{t). Then we have 

(6) nB(o)(e,r?) = i^- / Hi.'n.t)dt. 

Jo 

Example 2.3. Let us say that a simply connected bounded domain B is 
algebraic of degree d if its Cauchy transform hs = —z + 4:dzIlB (which is 
analytic in B) is actually a rational function of degree d — 1 (see \Gu\ lEV] ). 
In this case the same is true for B(t) for t > 0, and thus it may be shown, 
similarly to |Guj . that the boundary of B{t) is defined by the algebraic 
equation Q{z,z) = 0, where Q is a polynomial of degree 2d. The genus 
of this algebraic curve is thus < {2d — l){d — 1), and thus the number of 
components of B(t) is at most {2d — l){d — 1). 

2.4. Points of complete contraction. Let us say that a point of is a 
point of complete contraction if it belongs to B{t) for all t E [0, t*). Thus 
the points of complete contraction are the points of disappearance of the 
bubbles which "survive" until the time of complete contraction t*. The set 
of all points of complete contraction is the intersection rit<:t*B{t). In |EE1] 
we described the structure of this set. 

Theorem 2.4. (i) The set of points of complete contraction is finite. 

(a) The points of complete contraction are the global minimum points of 
the gravity potential. 

Remark 2.5. Note that if the boundary of a domain B is smooth, then by 
the Hopf's strong maximum principle, a global minimum of II^ cannot be 
attained on the boundary of B. 

2.5. Asymptotics of contraction. It is shown in [EVj that when a bubble 
contracts completely to a point, and the Hessian of the potential at that 
point is nondegenerate, then the boundary of the bubble has the asymptotic 
shape of an ellipse, whose half- axes are directed along the eigenvectors of 
the Hessian, and their lengths are inverse proportional to its eigenvalues. 
Here we would like to extend this result to the case when the Hessian may 
be degenerate. 

Namely, assume that the bubble B contracts at a point 0. In this case, 
IIb{0) is an isolated global minimim point of H^, [EEll lEV) . Let us assume 
that is a degenerate critical point, and the kernel of the Hessian of H^ at 
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zero is the x-axis, i.e. IIb{x, y) = + 0(|2;|^) near 0. For simplicity let us 
first assume that the bubble is symmetric with respect to the x-axis. Then 

where (3 is some positive number, and ... are monomials strictly inside the 
Newton polygon. Let us call n the degree of the critical point 0. 

Let B^:{t) be image of the bubble at the time t under the renormalization 
X — > cx, y — > c^^'^y, where c = c{t) is chosen in such a way that the 
diameter of B^,{t) is 2 (so c{t) behaves like {t* - t)"^/^" as t ^ t*). 

Define the polynomials 



Qniu) = J2 



Ml 

k=0 



U 



Theorem 2.6. The boundary of the domain tends to the curve 

y^ = (3\l-x^)Ql_,{x''). 
Proof. We write in the the form 

Ub = -liz - zf + /3Re(z2«/2n) + 
o 

where ... are the terms strictly inside the Newton polygon. 

The conformal map of the unit disk into the outside of B(t) which maps 
to oo has the form 

f{c) = A{c+c') + m, 

where (p is an odd holomorphic function in the unit disk, and one may assume 
that ^ > 0. Here (j) = c^t, A = A{t), and A ^ as t ^ t* . 

Let /* denote the function obtained from / by conjugating the coefficients 
of the Taylor series, and let hsiz) be the Cauchy transform of B: 

hB = z- AO.Ub = z- 2/3z2"-i + ©(z^"). 
Then by Richardson's theorem (see |R H lEVj ) . the function 

f*{C')-hB{f{0) 
extends holomorphically from the boundary of the unit disk to its interior, 
and vanishes at 0. Therefore, we obtain, for |(^| = 1: 

Im(/)(C) = -2pA^''-hm{C + C"^)+ + 0{A'^''), 

where the subscript + stands for the holomorphic part (i.e. the nonnegative 
degree terms of the Taylor series) . 
Now we claim that 

Im(C + C"')+ = 4"-iQ„_i(cos2 e) sine, 

where = e*^. This is proved easily by induction in n. Also, if |C| = 1 then 
C + is real. Thus, we get for |C| = 1: 

y = Im/(C) = -/3(2A)2«-iQ„„i(cos2 0)sin0 + O(A2n)^ 
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On the other hand, we have 

x = Re/(C) = cos + 0(^2). 

Therefore, upon rescahng x — > x/2A,y 7//(2A)-^"~^ we wih obtain the 
fohowing equations: 

y = -/3g„_i(cos2 0) sine + 0(A), 
X = cosd + 0{A), 
Thus when A goes to zero (i.e. for t ^ t*), we get the hmiting shape 

y^ = P\l-x^)Ql_,{x^), 
as desired. □ 

Now let us consider the general bubble which is not necessarily 

symmetric with respect to the x-axis. In this case, 

UB = \y' + ayx^ + ^x^^ + ...= 
2 In 

-^(^ - + Re(ia— - + /?— ) + 
8 n + 1 2n 

where a, f3 are real and P — na^ > (to insure that is an isolated minimum 

of Hb), and 

hB = z + 2iaz^ - 2/3z2"-i + 0{z^''). 
Applying a similar method to the one used in the symmetric case, we obtain 
the following result. 

Let B^{t) be image of the bubble at the time t under the map z — > z+iaz"^, 
followed by the renormalization x cx, y — > (?'"'~^y, where c = c{t) is 
chosen in such a way that the diameter of -B*(t) is 2. 

Theorem 2.7. The boundary of the domain B^{t) tends to the curve 
y^ = [P-na^f{l-x^)Ql_^{x^). 

The most interesting case for applications is n = 2. This case corresponds 
to contraction of (symmetric) domains that are on the boundary between 
those that break up and those that don't. 

Example 2.8. Assume that = ^+/?Re2;^/4. This potential corresponds 
to the contracting bubble whose conformal map from the unit disk to the 
outside region has the form 

where A = A{t) > is some function. (This map is found from the singu- 
larity correspondence, |E Vj ) . For small enough A, this function is univalent 
and defines a bubble. Contraction of the bubble corresponds to decreasing 
A to 0, leaving (3 fixed. Then the bubble contracts to the origin with the 
asymptotic shape given by the above theorem for n = 2. This domain is on 
the boundary between rupturing and non-rupturing domains. 
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Remark 2.9. This analysis of asymptotic shapes is similar to the analysis of 
the shapes of the necks of bubbles during break-off which is done in [LBW| . 



3. Points of partial contraction 

3.1. Definition and properties of points of partial contraction. Let 

B()(t) C B{t) be a connected component of the air domain (i.e. a single 
bubble), which exists on the interval of time (T/,rc); namely, we assume 
that Tf is the time of formation of the bubble Bq, and Tc is the time of its 
disappearance (contraction). A point contained in Bo{t) for all t € {Tf,Tc) 
will be called a point of partial contraction. In Fig.l, P is a point of com- 
plete contraction, and Q is a point of partial contraction. A point of either 
complete or partial contraction will be called a contraction point. To every 
contraction point there corresponds a time of contraction Tc. 

Theorem 3.1. (i) Every component Bq of the air domain which contracts 
without breakup contains a unique contraction point. 

(ii) A contraction point at a time t is an (isolated) global minimum point 
of the potential 'i^B{o)\B{t) ^ '^^'^ ^^ce versa. In particular, the number of 
contraction points is finite. 

Proof. By formula 

ns(o)\iJ(t) (C, v) = Kt- I $(^, r?, T)dT, 

Jo 

where Kt is a constant. Since <I>(^,?7, r) < 0, and ^{S,,r],T) = if and 
only if {(,,r]) is contained in the closure of B{t), we have ^B(o)\B{t) ^ ^t, 
and n5(o)\B(t) = Kt if and only if (^,77) G r\r<tB{T). Thus, if (^,r/) is a 
contraction point at a time t, then ^B{o)\B{t)iCi'n) = ^t, i-e. the potential 
achieves its minimal value at {(,,rj), and vice versa. 

On the other hand, the contraction points of the component Bq are con- 
tained in the domain B{0) \ B{t), where the potential 'n.B{o)\B{t) satisfies 
the Poisson equation AH = 1, i.e. is a real analytic function. Thus the set 
of contraction points of Bq is analytic (as it is a connected component of 
the set of solutions of the equation n^(o)\_B(t) = Kt). Also, it is compact 
and simply connected. This implies that it consists of one point (see [EEIJ ). 
This point is thus an isolated point of global minimum of the potential. On 
the other hand, it is clear that any global minimum point of ^B{o)\B{t) is a 
contraction point. The theorem is proved. □ 

3.2. Finding points of partial contraction. If the bubble breaks up 
into two parts, which contract without further breakup, then Theorem 13.11 
allows us to find the partial contraction point explicitly, provided that the 
derivative of the conformal map /o(C) from the unit disk to the complement 
of the initial domain -B(O) is a rational function. Indeed, let r be the moment 
of contraction of the bubble that contracts sooner, and zq be its contraction 
point. The domain B{t) is connected, so there exists a conformal map 
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friC) of the unit disk into the complement of B{t), which also has rational 
derivative. It can be found as described in [EVj . Assume that the map fr is 
known. Then by a direct computation one finds the potential n^(^) on the 
whole plane. Next, the point zq is found from the condition 

v(nB(o)-nBM)(zo) = o. 

Moreover, the value of the potential n^(Q) — n^j-^) at zq must coincide with 
the value of this potential inside the domain B[t) (this value is constant 
inside B{t) by Theorem I2.2p . This is a condition on the unknown time r of 
contraction. 

A similar method allows one to find the contraction point of the smaller 
bubble in the problem of simultaneous contraction of two circular bubbles, 
Ki = {z; l^l < R} and K2 = {z; — a| < r}, a, i?, r > 0, i? > r, a > i? + r 
(Fig.2). 

In this case, as was shown by P. P. Kufarev |Kuj . the conformal map of 
the unit disk to the complement of the domain B{t) at the time r of partial 
contraction has the form 

^-(^) = T3^+C' ^=2l^"" + ^ )' 

o 2r2 + i?2_^\ 

P = aa ^ 

2 y aa J 

where a > and is the middle root of the cubic equation 

(8) 2a^x^ - (2a2(i?2 _ qt/^^^ + a'^^x^ + {2r^ + - qt/^f = 0, 
and from the equations 

(9) VnB(o)\B(T)(2:o) = 0, '^B{a)\B(T){zo) ='^B{a)\B{T)\B(T) 

one finds the contraction point zq of the smaller bubble, and the time r of 
partial contraction. 

3.3. Asymptotics of partial contraction. As in the case of complete 
contraction, if the contraction point is a nondegenerate minimum of the po- 
tential, then at times close to the time r of partial contraction, the boundary 
of the bubble has an approximate shape of an ellipse, whose half-axes are 
directed along the eigenvectors of the Hessian of I[b{o)\b{t) the contrac- 
tion point, and their lengths are inverse proportional to its eigenvalues. This 
fact is proved analogously to the case of complete contraction ( |EVj ). 

It is interesting to study the rate of partial contraction, in the case when 
the bubbles contract at different times. Namely, assume we have two bubbles 
Bi and B2, and B2 contracts to a point P2 at a time t', at which Bi assumes 
the shape of a domain E. We assume that P2 is a nondegenerate minimum 
of the potential ^B{<d)\B{t')- Let us conformally map the outside of E onto 
the unit disk, so that 00 maps to 0, and the point P2 to a point h G (0, 1). 



Obviously, such b and the map are unique. Let ^ be the complex coordinate 
in the disk; then the potential $ at a time t < t' close to t' has the form 



* = -^l°slCl + |loglf3^1 + O(t'-0, 

where Q = Q{t) is the rate of contraction of the bubble B2; it is clear that 
Q{t) -^Oast^t'. 

Let T{t) be the image of the boundary of B2{t) in the disk, and let A{t) 
be the area of B2{t). Since the boundary of B2{t) for t close to t' is almost 
elliptic, the distance from the points of T to the point b is sandwiched 
between ciA{t)^/'^ and C2^(t)^/^, where ci,C2 are some constants. On the 
other hand, $ must vanish on T{t). This yields 

Q , 

— \ogA^q\ogb,t^t. 

Thus, Q{t) is equivalent to 2q iog"^(t) as t — > t' . 

This allows us to determine the asymptotic behavior of A{t) as t ^ t' . 
To do so, let us introduce the variable t = t' — t. Then asymptotially A 
behaves as the solution of the the differential equation 

dA ^ ^^ log6 
dr log A ' 

with the initial condition A = as r = 0. Solving this equation, we get the 
solution which is implicitly defined by the equation 

AlogA- A = 2qT log b. 

So we obtain the following result. 

Theorem 3.2. We have 

2qT\ogb 2qlogb 

A , Q , T^O. 

log r log r 

This result can be generalized to the case when the bubble B2 contracts 
to a degenerate minimum. Namely, assume that the contraction point is a 
critical point of degree n. Then, conducting a similar asymptotic analysis, 
using the results of subsection 2.5, we obtain the following theorem. 

Theorem 3.3. We have 

2nqT\ogb 2nq\ogb 

A^ — , Q^— , r^O. 

log r log r 

Thus we see that when t is close to t' then the contraction of the bubble 
B2 is logarithmically slow, and almost all air is extracted from Bi. 
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4. Sufficient conditions for breakup of symmetric bubbles 



Theorem 4.1. Assume that the gravity potential of a simply connected 
domain B{Q) symmetric with respect to a point P (respectively, a line £) 
achieves a global minimum at a point Q ^ P (respectively, Q ^ i). Then 
B{0) breaks up in the process of contraction. 

Proof. The symmetric point Q' ^ Q\s also the global minimum point for the 
gravity potential. By Theorem[3l the points Q and Q' are points of complete 
contraction of the domain 5(0). Therefore, -6(0) must break up. □ 

Corollary 4.2. (see [EElj ) Let B{0) = {{x,y);\x\ <b,y'^ < f{x)}, where 
f{x) is an even piecewise smooth function, nonnegative on [—b,b], such that 
/(6) = (Ftg.3). 

Then B{0) breaks up in the process of contraction if 

, , f^d(xf(xy/^) 
10 / \^ L / >vr/2. 

Jo + fix) 

Proof If ([lOD holds, then 



dx"^ 2 TT Jq x"^ + f{x 

so the origin is not a global minimum point of the potential. Hence the global 
minimum is attained at another point. Because of the central symmetry, the 
domain B{0) must break up. □ 



Theorem 4.3. (see [E Vj . problem 2 on p. 32) Let B{0) be a simply connected 
domain, symmetric with respect to the horizontal axis, and the function 
7r(x) := 115(0) (^)O) more than one local extremum on (— cxd,cxd). Then 
the domain B{0) breaks up in the process of contraction. 

Proof. Assume the contrary, i.e. that the domain does not break up. Let 
xi{t) < X2{t) be the intersection points of the boundary dB{t) with the 
horizontal axis (Fig. 4) 

By the monotonicity property of contraction, the function xi{t) is in- 
creasing, and the function X2{t) is decreasing on [0,t*), and limj^t* xi{t) = 
limf^i* X2{t) = xq. The rays [x2{t), +00) and (—00, a;i(t)] are flowlines of the 
flow, and the flow is directed from infinity, so the potential <l>(x, 0, t) increases 
from —00 to on the interval (— oo,xi(t)], equals zero on {xi{t) , X2{t)) , 
and decreases from to —00 on [x2(t), +c>o). Thus, if ii > £,2 ^ xq or 
6 < 6 < Xq, then for any t £ {0,t*) one has $(^i,0,i) < ^>(6,0,t). More- 
over, if this inequality turns into an equality for all t, then ^1 = (,2- Since 
n5(o(^,r/) = K — Jq ^{^,r],t)dt, these arguments imply that the function 
IIb(o){x,0) is strictly increasing on (xo,+oo), and strictly decreasing on 
(— oo,xo), i.e. its unique local extremum is a minimum at the point xq. A 
contradiction. □ 
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5. Regulated contraction 



5.1. Definition of regulated contraction. Consider contraction of a do- 
main which breaks up into two parts at a time r E [0,t*). After the time 
T, the process of contraction may be controlled, creating different pressures 
inside the two components of the air domain by regulating the amount of 
air which is pumped out of each component. This is a generalization of 
the problem from Section 2. In particular, it is interesting whether one can 
regulate contraction in such a way that both bubbles contract at the same 
time; in this case, as we will see, the contraction points have the convenient 
property that they are critical points of the potential of the initial bubble, 
and thus can be easily found. Let us consider this generalized problem in 
more detail. 

Consider the process of contraction of a system of two bubbles; we don't 
assume that they were obtained as a result of breakup of a single bubble. 
Assume that at a time t S [0, i*), the air domain B{t) consists of the com- 
ponents Bi{t) and B2{t), and the air is pumped from Bi{t) at the rate 
qi{t), and from B2{t) at a rate q2{t). This means that the velocity potential 
^>(x, y, t) is a solution of the boundary value problem 

^+5^ = 0' (^,y)ei?(tr; $|aB,(i) = $^W,i = l,2; 

^(x, v) = -SM±Ml iog(,) + 0(1), , ^ 00, 
zvr 

and the constants and $2(i) are chosen in such a way that 

(11) / —di = q^{t),i = 1,2, 



dn 



It is useful to extend ^> to the interior of B{t): ^ = $i(t) in Bi{t) for i = 1,2. 
Then $(x, y, t) is an everywhere continuous function. The velocity of motion 
of the boundaries dBi and dB2 is The motion can be considered up 
to the time of disappearance of one of the bubbles (we assume that there 
is no topological transformations of the first kind, i.e., formations of new 
bubbles). The contraction process in this situation will be called regulated 
contraction (as opposed to free contraction, i.e. with equal pressures in the 
bubbles). The vector- function (qi{t), q2{t)) will be called the strategy of air 
extraction. 

Theorem 5.1. The gravity potential inside every component of the contract- 
ing domain changes by a constant in the process of regulated contraction ( the 
constants may he different for different components). 



This theorem is proved analogously to part 1 of Theorem 12.21 

Theorem 5.2. Let he two smooth families of simply connected do- 

mains, s = {si,S2), Bf n B2 = 0. // the potential of the system {i?f,i3|} 
changes hy a constant inside Bf and inside under the change of s, and 
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the areas of Bf and i?| are independent of s, then Bf and i?| don't change 
under the change of s. 

Thus, a pair of simply connected domains is locally uniquely determined 
by the potential inside them, given up to additive constants, and by the 
areas of the domains. 

Proof. Let Vg : dBf U 5S| ^ M is the velocity of the boundary under the 
change of s. For any continuous function u{x,y), 



Thus, the conditions of the theorem imply that the function Vg is orthogonal 
on dBf U to the following functions: 

(i) the characteristic functions £,j{z) of the boundary comnponents dBj 
(because the areas of Bj are constant in s); and 

(ii) ri^{z) = Re{w — z)~^ , 77^(2;) = Im(u; — z)~^, where w £ Bf U B2 
(because of the fact that the gradient of the gravity potential is constant in 
s). These functions are dense in L'^{dBf U Si?!)- Hence, Vs is identically 
zero. The theorem is proved. □ 

Theorem 5.3. Let qs{t), t G {t,0), Q = {<li,Q2), s G [0,1], be a smooth 
family of strategies of extraction such that qs{t)dt = Q is independent of 
s. Let B'^{9) = Bf{6) U i?!!^) result of extraction according to the 

strategy Qs from the same initial domain B{t) = Bi{t)UB2{t). Then B'^{9) 
is independent of s. 

Proof. By Theorem 15. 11 the gravity potential B^{9) is independent of s (and 
equals the potential of B{t)) up to an additive constant (in each connected 
component). The areas of the domains Bi{9) and -BK^) are also independent 
of s and equal Si — Qi, S2 — Q2, respectively, where Sj are the areas of 
the components Bj{T), and Qj are the coordinates of the vector Q, i.e., 
the volumes of the air extracted from the first and the second bubble. By 
Theorem 15.21 Bj{9) do not change under the change of s, as desired. □ 

Thus, the result of contraction depends only on the total quantities of air 
extracted from the bubbles for a given period of time, and does not depend 
on other parameters of the strategy. In other words, the transformations 
in the space of domains defined by extraction of air from the first and the 
second bubble, respectively, commute with each other. This is an analog 
of Richardson's result |Ri] on the commutativity of injection operations at 
different points (see also |EVj ) . 

5.2. The phase rectangle and the accessibility region. Theorem 15.31 
shows that the domains which can be obtained from B{0) under regulated 
contaction can be visualized by points of the "phase" rectangle < X < Si, 
< Y < S2, where Si and 5*2 are areas of Bi{t) and -82(7"), respectively: to 
the point (X, Y) corresponds the domain which is obtained by extraction of 
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the volumes Si — X,S2 — Y from the first and second bubble, respectively. 
The strategy of extraction is depicted by a path inside the rectangle, which 
emanates from the corner (5*1,52) (Fig. 5). 

It is important to note that in the process of regulated contraction, pieces 
of the boundary may move towards the fluid region. If this happens for a 
particular strategy, the problem of regulated contraction is ill-posed in the 
vicinity of this strategy. Thus in general one has neither the monotonicity 
property nor the existence of a weak solution up to the time of complete 
contraction. In other words, the solution of the problem of regulated con- 
traction may not exist for some strategies of extraction. More specifically, 
in a generic situation solutions develop singularities in the following way: 
at some point of the boundary one of the two bubbles develops a semicu- 
bic cusp directed towards the fluid region (Fig. 6), and after this time the 
solution cannot be continued. 

The development of such a cusp was first discovered by P. Ya. Polubarinova- 
Kochina |Ko| for the problem of contraction of the boundary of the oil region. 
The fact that the singularity generically takes the shape of a semicubic cusp 
is related to the fact that while the solution exists, the boundary of the 
bubble is an analytic curve, and a semicubic cusp is the simplest singularity 
of a nonselfintersecting analytic curve. 

It follows from the above that not all paths in the phase rectangle corre- 
spond to actual solutions, but only those that lie in some region O, which 
is the set of all points of the rectangle that can be accessed by a strategy of 
extraction in which both bubbles exist all the way up to the last moment. 
We will call the accessibility region. The boundary of the accessibility re- 
gion consists (in the generic situation) of parts of the boundary of the phase 
rectangle, and curves, whose points correspond to pairs of bubbles, one of 
which has a cusp (Fig. 7) 

The path 7 corresponding to free contraction divides the accessibility 
region into two parts. Namely, if the pressure in the first bubble is kept 
higher than in the second bubble, then the corresponding path lies below 7, 
and if it is kept lower then the path lies above 7. Similarly, one can construct 
the trajectory 7p of free contraction, starting from a domain corresponding 
to any point P G 17. The accessibility region Vt is foliated by such trajectories 
(Fig. 8), which implies that is contractible (i.e., simply connected). 

If the initial domain consists of two symmetric bubbles, then the accessi- 
bility region is symmetric with respect to the diagonal Y = X oi the phase 
square, and the path 7 of the free contraction is this diagonal. 

Remark 5.4. Theorem 13.21 implies that the trajectories 7p of free contrac- 
tion which do not end in the origin are tangent to the boundary of the 
phase rectangle at the endpoint, and the tangency is of the type Y — a = 
cX\og{\/X) if the endpoint is (0, a), and X — a = cyiog(l/y) if the end- 
point is (a,0) (for some c > 0). 
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6. Synchronizing strategies of extraction 

Let us say that a strategy q{t) is synchronizing for the system of bubbles 
Bi, B2, if the extraction according to this strategy leads to simultaneous 
contraction of the bubbles Bi and B2 to a point. The path in the accessibility 
region which corresponds to a synchronizing strategy ends in the origin. 
Obviously, a synchronizing strategy exists iff the accessibility region contains 
the origin. A domain B(0) which breaks up under contraction into two 
bubbles which have this property, will be called synchronizable. 

Under extraction of air according to synchronizing strategy q(t), the bub- 
bles Bi and B2 simultaneously contract to points Pi and P2. This means 
that the points Pi and P2 are limiting positions of the boundaries dBi{t) 
and dB2{t) when the time t tends to the time t* of contraction. These points 
are easily found from the shape of the initial domain. 

Theorem 6.1. The contraction points Pi and P2 are critical points of the 
potential ^b{o) (if they both belong to B{0)). 

Proof. 

Lemma 6.2. Let G be a bounded domain of area S. Then 
(12) \VUG{x,y)\<y^,x,y£R. 

Proof. Let z = x + iy,w = u + iv. Let K be the disk of area S centered at 
{x,y). Then we have 

,„„ , 1 , f {z — w)dudv , 1 f dudv 

|vnc(.,rt| = -I ^^YJaV^\- 

2'kJk\z-w\ 27r7o Jo 
as required. □ 

Now we prove the theorem. Let t„ be a sequence of times which tends 
from below to t* . Let P^\ P2"'^ be sequences of points in Bi{tn) and B2{tn), 
which converge to Pi and P2 as n — 00. By Theorem 15.11 we have 

vns(o)(P,("))=vns(,„)(p,(")). 

By Lemma [6^ \'VIlB[t„)(.Pj^^)\ — > as n — > 00, as the area of B{tn) tends to 
for n ^ 00. Hence, Vn^(o)(^j-"^) ^ 0, n ^ 00, j = 1, 2. Since the gravity 
potential is a C^-function, this implies that Vn5(o)(-Pj) = Oj ^ desired. □ 

Clearly, one of the contraction points Pi,P2 coincides with the point P 
of complete contraction in the sense of Section 1 (namely, the contraction 
point for the bubble that contracts later under free contraction). This point, 
as we mentioned, is the global minimum point of the gravity potential. Re- 
garding the second point, it is shown below that it is either a local minimum 
point or a degenerate critical point, and the degenerate critical point arises 
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for initial domains which he on the boundary between synchronizable and 
nonsynchronizable domains in the space of domains. 

Let us call an initial domain -B(O) strictly synchronizable if the corre- 
sponding accessibility region contains a sector 

{{X, Y)\X >0,Y> 0, +Y'^ < e} 

for sufficiently small e. Clearly, strict synchronizability is an open condition, 
i.e. this is a property which is stable under small deformations. Non-strictly 
synchronizable domains form the boundary between synchronizable and non- 
synchronizable domains. To illustrate this, we show in Fig.9,a,b,c what the 
accessibility region looks like for a strictly synchronizable, non-strictly syn- 
chronizable, and nonsynchronizable domain. 

Theorem 6.3. // a domain B{0) is strictly synchronizable, then the con- 
traction points Pi, P2 for a synchronizing strategy are local minima of the 
potential B{0) (if they lie in B(0)). 

Proof. If the domain B{0) is strictly synchronizable, then one of the trajec- 
tories of free contraction ends in the origin (Fig. 9). 

Therefore, there exists a synchronizing strategy which corresponds to free 
contraction on the interval {t* — e, t*) for some e > 0. Then by Theorem l3.lt 
the points Pi,P2 are points of global minimum of the potential IiB{t*-£)- 
Since inside Bi{t* —e),i = 1,2, the potential Ilt*-e coincides with n^^g) up 
to constants, we see that Pi and P2 are points of local minimum of IIb(o) ■ ^ 

Remark 6.4. If the point Pi or P2 is outside the domain B(0) (a priori, such 
a situation cannot be ruled out because of the failure of the monotonicity 
property), then it is a local minimum point of the analytic continuation of 
the potential of of the domain B{0) from its inside to its outside along the 
track of the corresponding bubble. 

7. ASYMPTOTICS OF REGULATED CONTRACTION UNDER A 
SYNCHRONIZING STRATEGY 

Theorem 7.1. Let B[t) = {Bi{t), B2{t)) be the evolution of the air domain 
after its breakup, under extraction of air according to a synchronizing strat- 
egy q{t), Pi,P2 G B{0) be the contraction points, and Ai,A2 the Hessian 
matrices of the potential IiB{o) these points. Let Bj{t) be the domains 
obtained from Bj{t) by dilation Xj{t) times, where Xj{t) is chosen in such 

a way that Bj{t) has a fixed diameter d. Then if the matrices Ai,A2 are 

nondegenerate, then the boundaries dBi{t) tend to ellipses, whose axes are 
directed along eigenvectors of Ai,A2, and the lengths of the half-axes are 
inverse proportional to the eigenvalues of these matrices. 

Proof. We prove the theorem for Bi; the proof for B2 is the same. Assume 
that the point Pi is the origin. Obviously, for t close to t* the potential 
11--— has a nondegenerate local minimum at some point a{t) inside Bi{t). 
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Let E{t) = Bi{t) — a{t) be the domain obtained by translating Bi{t) by 
the vector —a{t). It can be seen, along the lines of [EVj . that the boundary 
dE[t) converges to some curve V. 

Let A£,[z) be the matrix of second derivatives of the gravity potential of 
a domain D. It is easy to see that AE(t){z) — > ^s(o)(0) = Ai, t ^ t* , M 
z S E{t) for t close to t* . This implies that the potential of the domain E 
bounded by the curve F is a quadratic function, whose Hessian matrix is 
Ai, and which has a minimum at zero, i.e. it is {Aiz,z) + C. By Sakai's 
theorem, this implies (see e.g. |EVj ) that E is an ellipse, whose axes are 
directed along the eigenvectors of Ai , and the lengths of half-axes are inverse 
proportional to its eigenvalues. The center of the ellipse is situated at the 
origin. The theorem is proved. □ 

8. Properties of the potential of a non-strictly 
synchronizable domain 

8.1. Contraction at degenerate critical points. 

Theorem 8.1. Let B{0) be a synchronizable domain, Pi and P2 he points 
of its contraction under a synchronizing strategy, and Ai,A2 be matrices 
of second derivatives of the potential at the points Pi,P2. In this case, if 
the matrices Ai,A2 are nondegenerate, then the domain B{0) is strictly 
synchronizable. 

Proof. Assume the contrary, i.e. that the domain is not strictly synchroniz- 
able. Consider the synchronizing strategy, which corresponds to a path 7 in 
the accessibility region (Fig. 10). 

We may assume that the path 7 goes along the boundary of 17. In this 
case, the boundary of one of the bubbles has a persistent singularity in the 
process of contraction. But by Theorem l7.ll the limiting shape of this bubble 
is an ellipse, which does not have a singularity. This is a contradiction. □ 

Corollary 8.2. // a domain is synchronizable, but not strictly synchroniz- 
able, then the matrix A2 of second derivatives of the potential at the point 
P2 is degenerate (here Pi is the point of complete contraction, and P2 is the 
second contraction point). 

In the generic situation this means that the potential has at P2 a critical 
point of type saddle-node, i.e. in some Cartesian coordinates 

where ... stand for the terms inside the Newton polygon. 

8.2. Asymptotics of contraction. It is easy to show that in the situation 
Corollary [121 the limiting curve F is a line interval (slit) of length d (traveled 
forward and backward), in the direction of the kernel of the matrix A2. 
However, it is interesting to study a finer asymptotics of contraction in 
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this situation. Consider the generic case when the contraction point of the 
singular bubble B2{t) is 0, and the potential at this point is as above: 

ni?(o)(^) g — + — 3 — + - 

Let B^:{t) be image of the bubble i?2(i) at the time t under the rcnormal- 
ization x —>■ cx, y — > c^y, where c = c{t) is chosen in such a way that the 
diameter of B^,{t) is 2. 



Theorem 8.3. (i) The boundary of B^{t) tends to the curve 

y^ = (3^(x + l-) (x-^ 



2J \ 2 

(a) The bubble Bi{t) contracts at some point a < of the real axis. 
Thus the cusp of B2 is always directed precisely towards Bi at the time of 
contraction. 

Proof. The method of proof is the same as the one we used for the asymptotic 
analysis of a bubble contracting to a degenerate minimum. Namely, we have 

hB,it)iz) = K{t) + z-2Pz^ + ..., 

where ... are negligible terms as t ^ t*. (Note that hB2{t){z) depends on t 
because of the presence of the second bubble Bi{t)). Thus, the conformal 
map of the unit disk into the outside of Bi (t) which maps to 00 looks like 

/(C) = ^(C + C') + Be + cc + D + ... 

Since /*(C~^) — ^d(/(C)) is holomorphic and vanishes at infinity, we get 

B = -2pA^,C = -AI3AD, 

modulo negligible terms. Also, we have a cusp in the "saddle" direction of 
the singularity, i.e. in our case in the negative direction. Thus, //(—I) = 0, 
which gives 2B — C = modulo negligible terms. This means that modulo 
negligible terms we have 

D = A, C = -A(3A^. 



This implies that 



Re/ = 2A{cose + ]^) + 0{A^) 



and 

y = -2l3A^{s\n2e + 2sine) + 0{A^). 

Thus after the rescaling x —>■ x/2A, y = y/^A^ and sending A to zero we 
get the limiting curve 
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and part (i) of the theorem follows. We also get K = QI3A^ + 0{A^) > 
for small A, which implies that the second bubble disappears at some point 



We see from the proof of this theorem that when A goes to zero, the 
area of Bi goes down as ciA^ and the area of B2 as 02-^ ■ This shows 
that the trajectory corresponding to our strategy (i.e. the upper boundary 
of the accessibility region) behaves near the origin as a semicubic parabola 



8.3. A sufficient condition of breakup. Corollary 18.21 gives us the fol- 
lowing sufficient condition of breakup. 

Corollary 8.4. Suppose , s £ [0, 1], is a smooth family of simply con- 
nected domains, such that B^ is symmetric (centrally or axially), and the 
potential of B^ admits a global minimum point in B^ which is not fixed by 
the symmetry. Assume that the analytic continuation of the potential of B" 
for any s does not have degenerate critical points. Assume also that it is 
known that B" undergo at most one topological transformation under con- 
traction for all s. Then the domains B^ break up under contraction for all 
s. 

Proof. By Theorem 14.11 B^ is synchronizable (and the synchronizing strat- 
egy is free extraction). If some B^ does not break up, then it is not synchro- 
nizable, so for some a < s B^ is synchronizable, but not strictly synchroniz- 
able. By Corollarv 18.21 the analytic continuation of XI^ct has a degenerate 
critical point. Contradiction. □ 

9. The boundary between rupturing and non-rupturing domains 

Having multiple local minima of the gravity potential is not a necessary 
condition for the breakup of a bubble. Indeed, it is obvious that a domain 
in Fig. 11 breaks up in the process of contraction, although its potential has 
a unique critical point (so this domain is not synchronizable). 

In this connection, it is interesting to study domains which lie on the 
boundary between rupturing and non-rupturing domains. 

Theorem 9.1. Suppose that B is a generic domain, which is on the bound- 
ary between rupturing and non-rupturing domains. Then B does not break 
up under contraction, but at some time t the domain B{t) has a cusp on 
the boundary, locally equivalent to y = x^/"^ . This singularity disappears 
immediately under further contraction. 

In other words, on the boundary between the sets of rupturing and non- 
rupturing domains in the space of (C'^-smooth) domains, a dense open set 
is formed by domains that have the property stated in the theorem. 

Proof, (sketch) Consider a smooth family of simply connected bounded do- 
main of generic position, B'^ , s G [0, 1]. Assume that for s < a the domain 

18 




□ 



Y = cX3/2 (Fig.9(b)). 



does not break up, while for s > a it does. Let t(s) be the time of 
breakup of the domain for s > a. Let r = lims^o-r(s). One can show 
that in the situation of general position this limit exists and is not equal to 
zero. Consider the family of curves r(s), s G (a, 1], which are obtained from 
B^ by contraction during the time t(s) (where t(o") := r)). 

The curves r(s) for s > a have a simple self-tangency at some point. It 
is easy to see that typical degenerations of such curves into simple closed 
curves have the structure described in the theorem: for s = a at the point 
of disappearance of the loop there forms a cusp of degree 5/2. Namely, a 
typical such family is 

= x'^{x + e),e = e{s), 
where e{a) = 0, and e(s) > for s > a. □ 

Remark 9.2. The appearance of the singularity in the process of contrac- 
tion is at first sight a strange phenomenon, as the contraction problem has 
good properties of correctness and stability. Nevertheless, solutions with 
such type of instantaneous singularities do exist. They were first discovered 
by Howison in |Ho| . More precisely, he showed that in the process of con- 
traction there can appear instantaneous cusps of degrees (4n -|- l)/2, n > 1, 
while cusps of degree (4n — l)/2 cannot appear. For n = 1, these are the 
instantaneous cusps of degree 5/2 that we have just considered. Solutions 
with these properties form a subset of codimension 1 in the space of all 
solutions. 

If we restrict ourselves to polynomial domains of degree < n, then the 
boundary between rupturing and non-rupturing domains, as follows from 
Theorem 19.11 is a piece of an algebraic surface in the space of coefficients. 
For small degrees n this equation is not difficult to write down. If two 
domains can be connected by a curve not intersecting this surface, then 
either they both break up under contraction, or they both don't. 

10. Corrections to [EETI IEV] 
We use the opportunity to correct some errors in [EE 11 lEVj . 

2 2 

1. The formula for the gravity potential of the ellipse ^ + p- = 1 given in 
[EETj lEV] ([EEl],p.517, and [M], (4.14)) is incorrect. The correct formula 
is 

2 a + b a + 

(i.e., a and b need to be switched). The same correction needs to be made in 
formulas (24), (25) of jEElj . and the two sentences after formula (25). Hence 
the correct condition of division of the bubble for a = 26 is c < ^/Sb. The 
same corrections should be made in [EVj in the example on p. 52, and in the 
answer to problem 4 on p. 68. 

2. As a result of 1, Theorem 4.10 in [ EV J should say that the lengths of the 
half-axes of the limiting ellipse are inverse proportional to the eigenvalues 
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of the Hessian matrix. The same correction is to be made on p. 535 of |EElj 
(the power —2 of the Hessian matrix should be replaced with +2). 

3. The left hand side of formula (19) should read 2^E^(0,0) (the factor 
of 2 is missing). The same correction should be made in the first formula 
on p. 52 of [EV ], The right hand side of the formula in Theorem 6.4 in |EVj 
should be tt, not 7r/2. In the first formula on p. 52 of |EVj . the factor 2/7r 
should be l/vr. 
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